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Abstract. Given a 2-dimensional no-pair Weyl operator Wz with a point 
nucleus of charge Z, we show that a homogeneous magnetic field does not 
lower the critical charge beyond which it collapses. 

1. Introduction 

Perfect graphene is modeled in continuum one-particle approximation by a two- 
dimensional free Weyl operator (massless Dirac operator). Non-perfect graphene 
' (— I '' has additional potentials; a particular case of importance is the presence of an 

impurity of Coulomb type (see the review of Castro Neto et al. [3]). As opposed 
to non-relativistic mechanics, in relativistic mechanics both kinetic energy and the 
Coulomb potential energy have the same linear scaling for large momenta which 
implies the existence of a critical coupling constant. This explains the interest 
in the subject in the physics literature, see, e.g., Pereira et al. )15| and Shyvtov 
et al. [17] . The critical coupling constant as occurring in these papers can be 
mathematically thought of as that coupling constant were a natural definition of 
self-adjointness ceases to exist. In addition to the electric impurity potential it is 
often also important to study the systems with an additional homogeneous magnetic 
field perpendicular to the graphene sheet. Of course, the question arises to what 
extend the presence of the magnetic field changes the critical coupling constant. 

If one is interested in multi-particle effects it is essential to have a well defined 
multi-particle Hamiltonian (see [5] and the references therein). Because of the 
Weyl operator's unboundedness from above and below, a naive addition of the 
one-particle operators acting on the various particles plus their interactions - as 
would be natural in non-relativistic quantum mechanics - does not give meaningful 
Hamiltonians (Brown and Ravenhall [2]). This problem can - on a physical level 
- be overcome by a quantum field theoretical treatment. Approximately, one can 
use the no-pair Hamiltonians initially introduced by Brown and Ravenhall [2| and 
further developed by Sucher jl8[ . Because a non-perturbative analytic treatment 
of quantum electrodynamics is not available, we will concentrate on the second 
alternative. 

A description of - one-particle ~ no-pair operators in an nutshell is as follows: 
the state space on which the no-pair Hamiltonians are defined depend on a Dirac 
sea in a similar way as the Fock representation of the electron-positron field de- 
pends on the initial splitting of the Hilbert space into electron and positron space 
(see Thaller [20l Section 10.1.1]). The Dirac sea is defined through an orthogonal 
projection (1 — A) in the state space of the Weyl operator Wz, i-e., the Hilbert 
space L^(C^,M^). The projection A is assumed fixed. The physically allowed states 
of Dirac particles will be those which are orthogonal to the sea, i.e., they are eigen- 
states of A. Metaphorically speaking the physical states are the vapor above the 
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Dirac sea. The no-pair approximation will then be the Weyl operator projected 
onto the states of fixed particle number iV - in our case A'^ = 1 - in the vapor, i.e., 
:= AWzA. 

It is reasonable to expect that expectations of the no-pair Hamiltonians Bz are 
bounded from below if Z is small (close to zero) and is unbounded from below for 
large Z. The critical coupling Zc constant is the value of the coupling constant 
where this change of behavior occurs. A priori Zc can be expected to depend on 
the choice of the Dirac sea. A particular simple choice is to take the Dirac sea as 
the one defined by the Weyl operator with the external homogeneous field. It is 
exactly this operator which we will be interested in. Our goal is to show that Zc 
does not depend on the presence of a homogeneous magnetic field. 

Although the rest of the paper is on the one-particle level, our interest in the one- 
particle stability stems from the multi-particle stability: the multi-particle energy 
of the no-pair Hamiltonian is bounded from below if and only if the corresponding 
one-particle Hamiltonian is bounded from below (see ). 

Our contribution is organized as follows: to escape the inconclusiveness of heuris- 
tic arguments, we give a precise mathematical formulation of the problem, collect 
some well known relevant facts and state our result (Section [2]). To prepare for the 
proof we study the partial wave analysis of underlying energy form in Section |3l In 
Section |4] we give the actual proof of our claim. The appendices contain auxiliary 
material which we collect for the convenience of the reader. 



2. Notation, Formulation of the Problem, and Main Result 

The Weyl operator (massless Dirac operator) of a particle of charge — e in two 
dimensions in a magnetic field diA2 — ^2^1 with vector potential 21 := {Ai,A2) 
and an electric potential ip is given by 

(1) Wa,<p := v(7-{p + ^21) - eip 

where v, c, and e are positive constants. Depending on the application, v could be, 
e.g., the velocity of light or the Fermi velocity in graphene, c is the velocity of light, 
and cr = (cri,cr2) are the first two Pauli matrices, i.e.. 

In this paper we are mainly interested in the case of a homogeneous magnetic field 
of strength B > orthogonal to the a;i-X2-plane, i.e., 2l(a;) = -j(— 2:2, xi), and 
an electric field generated by a nucleus of atomic number Z, i.e., ip{x) — Ze/\x\. 
(Note that B < corresponds just to a reflection of the coordinates x — >■ —x. For 
B = see Remark [5] of Section [21) This operator is to be self-adjointly realized 
in L^(M^,C^). Following Brown and Ravenhall [2] - see also Sucher [TSl [TH] - we 
will project these operators to the orthogonal space of a Dirac sea. More precisely, 
we are interested in the quadratic form of W<^^^ restricted to the positive spectral 
subspace 

^ := {V; e L^(M^C2) I V^ = AV'} 

with A := X(o,oo)(W2i,o)- 

By dilation W^ai.^ unitarily equivalent to a/ eB/ {2ch)vWc2z/(hv)- Thus it suf- 
fices to study Wz := cr ■ {p + {—X2,xi)) — j|j assuming that e = v ^ h = 1. We 
use complex notation z :— xi + 1x2, and, correspondingly d '■= ^{di + 1^2) and 
g l{di - 1^2) and introduce d := -2i(a+z/2) and d* := 2i{-d + z/2). This 
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allows us to write more compactly 



(2) 



Wz 



d 
d* 



i) 



Z 



We define qo as the linear span of the functions (f)m.n defined in (j36p of the 
Appendix. We also define Qq as the linear span of the spinors ipm,n defined in (|44|) 
of the Appendix. 

Theorem 1. The quadratic form {'4>,Wz'4') positive on Qq and extends to a 
closed form £ on Q which is bounded from below, if 



For Z > Zc the form is unbounded from below. 
We remark: 

(1) Physically £[il)] is the energy of an electron in the state ip on top of the 
Fermi sea defined by i^"*" in the quantum dot defined by the homogeneous 
magnetic field and an interstitial atom with charge Z. 

(2) li Z < Zc, then the form £ defines - according to Friedrichs [SI Satz 3] - 
a unique positive self-adjoint operator whose form domain includes qo and 
extends AWz- It is called the no-pair Hamiltonian of one electron in the 
quantum dot. 

(3) For scalar type magnetic operators, like Schrodinger operators (p — 21)^ -I- V 
or Chandrasekhar operators |p— 2t|-|-y, it is known that 2t does not lower the 
ground state energy because of the diamagnetic inequality. For operators 
involving spin in an essential way like the Pauli operator this is known to 
be false. Although, in our case, we cannot expect the energy to increase 
when 2t is turned on, our result shows, that the energy is not lowered 
dramatically, i.e., the critical coupling constant is not lowered. Thus, the 
boundedness result can be interpreted as a weak form of the diamagnetic 
inequality. 

(4) The result for Z > Zc means physically that the electron is pulled into 
the nucleus of the interstitial atom as the trial function of the proof will 
indicate. 

(5) The critical coupling constant in the three dimensional non-magnetic case 
with arbitrary non- negative mass was found by Evans et al [5] . Tix |211 122j 
sharpened the result to strict positivity with a lower bound linear in the 
mass. 

(6) The critical coupling constant in the 2-dimensional non-magnetic case was 
investigated by Bouzouina [I]. An error in the constant he obtained was 
corrected by Walter [53] . 

(7) The 3-dimensional magnetic case - for a rather big class of magnetic fields 
- was treated by Matte and Stockmeyer [13] . They showed that the critical 
constant is not lowered by an intricate resolvent method. The generality 
of their result is paid for by the absence of an explicit lower bound on the 
energy. The bonus of our direct approach based on Lieb and Yau's [13] 
strategy in the variant found in [5] - compared to transfering the methods 
of [2] - is our result on the positivity of the energy. 

(8) The numerical value of the critical coupling constant is Zc ~ 0.3780 which 
is compared with the expected critical coupling constant Zc of the existence 
of a distinguished self-adjoint extension of the non-magnetic Weyl operator 
Wz- Pereira et al [T5] and Shytov et al [T7] suggest in physical language 
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and using physical arguments that — 1/2. Recently Warmt [24, Satz 
2.2.6] showed that this is indeed the case. 



3. The positive spectral subspace and partial wave analysis 

The fact, that we are dealing with spinors in the positive spectral subspace of 
Wo allows us to reduce the problem to unrestricted scalar wave functions (see [51 
Section 1.1] for the three dimensional case). 

Lemma 1. The map 

(4) " 

is unitary. 

Furthermore, its restriction to qo is a unitary map from qo to Qo with the asso- 
ciated scalar products. 

Proof. First, we remark that $ maps indeed to Sj. This holds, since Sj is the closure 
of Qo in the L^-norm. 

To show that $ is surjective, assume -0 = {u, u)* e and orthogonal to 

{{w,d*\d*\~'wY\w e L2(M2)}. 

This implies 

{u,w) + {\d*\^^dv,w) = {u+ \d*\^^dv,w) = 

for ah w e L^{R'^), i.e., u = -\d*\-^dv. 
Next we remark that 



71 G No, m G Z are eigenvectors with negative eigenvalue, namely — 2^n + m+ + 1. 
Thus 

-\d*\-^dv\ ( 







for all n and m which implies {\d*\~^dv,4'rn,n) — 0, i.e., dv = 0. Therefore, -0 = 
(0, w)*. Such vectors are in the kernel of Wo, i-e., orthogonal to the positive spectral 
space, so that in the end = is the only vector in the positive spectral space 
which is orthogonal to <^{L'^{W' ,€?)). 

The identity {u, v) i2(R2-) = (<&m, for all u,v £ L'^iK) is immediate, as is the 
unitarity of the restriction. □ 

Using Lemma [1] we define the operator wz '■= ^*Wz^ on qo. The associated 
quadratic form on qo is 

(5) (it, wzu) := {u, $* Wz^m) = (u, \d* \u) - Z{u, Vu) 
with 

(6) V = \[y^+\d*\-'d^^d*\d*\-' 

Corollary 1. The operators AWz on Oo 0''nd Wz on qo are unitarily equivalent by 
Lemma [II In particular both operators and also the associated forms have all the 
same lower maximal bound. 
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Next we calculate the matrix elements of wz in the orthonormal basis given by 
the eigenfunctions (f>m,n of wq . First of all, we remark that this matrix is diagonal in 
the angular momentum quantum number m. We get for the matrix associated 
with Wq the following matrix elements 



(7) := {(j),n,n,Wo(f>m,n')Sn,n' = 2y/n + m+ + 1(5„,„' 

which is immediate from the eigenvalues equation (j45p : for the first summand of 
the potential V (see ([51)) we get the matrix w™'" with matrix elements 

7rV(" + l)|™|("' + l)|m| ^ (n-fc + i)i(n'-fc + i)i 

which is obtained by explicit calculation using the generating function of the gen- 
eralized Laguerre polynomials [12. Formula 22.9.15] and their recursion relations 
[T^ Formula 6.1.15]. (For convenience we use Pochhammer's notation {z)a 
T{z + a)/T{z) [see also (011)] ■) Eventually, the second summand of the potential V 
yields the matrix ti'"'^ with matrix elements 

1 ( m+l-O • r ^ n 

(9) t^rn' = (rf*M*rVm,n, J^d* \d* \-' <l>m.n') " " 



m+1,0 



if TO < 0. 



This can be obtained from ([5]) by observing that 

(10) d*Mrv™.4 l!"^-"' 

[-l(l>m+l,n+l if TO < 0. 

Thus, the quadratic form £„i of the matrix {4'm,mWz4'm,n'), for fixed angular 
momentum m G Z, on /o(P^o) ^ the subscript denotes sequences of compact support 
- is given as 

oo 

(11) = E ^ ['Zn' - f ««' + 

n,n'— 

As mentioned in Appendix El {4'm,n,wz(t>m\n') = (0m,n, w'z<^m,n')<^m,m', i-e., 
both, potential and kinetic energy, are diagonal in to. Thus, 

(12) {u,w,u)^Y.^^r.[a^^ 

where we write a™ := {(f>m,n,u) for the generalized Fourier coefficients for u S qo 
and where we collect those coefficients with the same angular momentum quantum 
number to and write 

(13) a'" = (a^,ar,---)- 
Obviously, {a'^)neNo. mez G ;o(N x Z). 

Lemma 2. The following facts for the matrix elements v^'^, of the Coulomb po- 
tential l/\z\ hold: 

• For all TO G Z and n, n' G No 

/-. ^\ n ^ m,0 |ni|,0 

1 14) < V , = V , . 



• For TO, 71, n' G No 



/I cr\ "^:0 \ m+1,0 

(15) V^ n' > 
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Proof. The first claim - including the remarkable positivity of all matrix elements 
- is immediate from the explicit expression ([5]). 

The second claim, i.e., monotony of the matrix elements in m, follows again from 
©, if 

r(fc + TO + 3/2) ^ r(fc + m + l/2) 



y^T{n + m + 2)T{n' + m + 1) ^r(n + m + l)T{n' + m + 1) 

for k < n,n'. This is immediate from the functional equation of the Gamma 
function. □ 

Lemma 3. We have 

(16) < {(f)m,n,V(j)„i^n') < (0O,n , 1^0O,«' ) 

for n, n' G Nq and m ^"L. 
Proof By ^ 



for all TO > 0. 

For negative to 



m,0 , m,l ^ 0,0 , 0,1 
V ' I + V , <■ V ' I + V I 



m,0 , m,l ^ -1.0 , -1,1 
V , + V , < V , + V , 
n.n' ' n.n' — 7i.n' ' n.n' 



where we use all three claims of Lemma [5J Thus, it sufhces to show that 

/I '7\ -1,0 I -1,1 ^ 0,0 , 0,1 

By © we have = and = ^'°'+i,„'+i- Thus ^ is equivalent to 

/io\ 0,0 ^ 0,0 

(18) ^^„+l,n' + l < 'f'n.n'- 

This is shown by induction, first in rt' and then in n. □ 
Corollary 2. for a G ^0(^0) 

we have 

inf{£:„[a]|a G ;g(No)} > M{£o[a]\a G ll(No)} 

and 

iui{{u,W:,u)\u G qo} = inf £o[^o(I^o)]• 
Proo/. Since the kinetic energy ^„ ^™n|an| is obviously invariant under the substi- 
tution a — >■ |a| and since the potential energy 

m , , m , 1 

7; ' -t- 77 

y "n,n' ' n.n' 

~^ / ^ On 2 

n,n' 

decreases by the same substitution because of the positivity of the potential matrix 
elements (Lemma [21 Formula (jl4p '). it suffices to take the infimum over positive 
sequences a G ^o(No) only. Thus, the desired inequalities follow from the corre- 
sponding inequalities of the matrix elements (jl6|) . □ 

4. Proof of the theorem 

Proof. (Theorem [T]) By Corollary [T] it is enough to study wz- By Corollary [2l this 
is equivalent to show lower boundedness of the quadratic form Sq on non-negative 
sequences a G ^^(No). 

At this point we embark on a strategy which goes back to Abel - at least - and 
which has been introduced in relativistic quantum mechanics by Lieb and Yau [13)1: 
it basically consists of estimating a non-diagonal operator by a diagonal one using 
the Schwarz inequality suitably. We will apply it to the two potential matrices 

and with matrix elements t'^'^/ and v^\, (jll[) . (For the matrix elements 
we will, from now on, suppress the reference to m = as well and write simply 
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wj^^/, (T e {0, 1}.) Given any sequence {gn)nefia with positive entries this strategy 
suggests estimating as follows: 



OO 2 

(19) ^a,v'^a)<Y,^Yl 

where we use matrix notation on the left and that v'' is symmetric. (Note that we 
suppress an index a with g although g can - and will - depend on a.) 
We start with the case a = and obtain 



(20) '"",tto " " to in-k+\).{n'-k + \). 

(fc + l)_. - 1 

n=0 fc=0 ^ 2 n'=0 ^ ' 2 

using and substituting n' — > n' + A;. We pick for a — 

(21) g„ ^ 



This allows to explicitly do the summation in n' and k which gives 

(00^ ( , ^ r(i)4 - r(n + |) , 

(22) („,,,) <__)^_^_^a„. 

n— ^ 4^ 

We apply Gautschi's inequality ([501) for G N and we get 

(23) £|i±|<^<^ 
which is also true for n = by inspection. Thus, 

(24) (a,«"a)<^^V^«^ 

ri=0 

It remains to treat the case a — \. We use again (fTO)) and obtain 

min{«V} (fc + l)i 



(a,v a) =— > — , > 
^ ^/(n + l)(n' + l) 



(25) „,,7loV(- + l)(-' + l) t^o (--fc+^)i(n'-fc+^)i 

^ h 9nV^ fro - ^ + fro ^n' + k + l{n' + i) . 
substituting n' — > n' + /c. In this case we pick 



(26) 5„ := 



which again allows for explicit summation in n' and k yielding 

32^2 ^ r(n + |) 2 



(27) («^-M<^E 



^4'' rt=o ^ 4 

For n > 2 we have 



(28) 



(.-i)-V2<Vm 
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By Gautschi's inequality (15(11) the left hand side majorizes r{n+ j)/T{n+ j). Thus, 
by the Gamma function's functional equation we get 

r(n + | ) 



(29) K < 



However, this inequality is also true for n = and n = 1 by inspection. Thus, 

QO_2 °o 

(30) ^a,v'a)<-j-J2V^al. 

''^i> n=0 

Putting all together we have 

oo 

(31) £o[a] > J2 2%A^n(l - Z/Zc)al > 



for Z < Zc- Note that the first inequality in ([3T|) is indeed strict unless a = 
because of ([^5]) . This shows the positivity of the form and therefore the first part 
of the theorem. 
This shows that 

(32) {u,v)q:^ ^ {{u,w^v) + {u,v)) 

nGNo, inGZ 

is a scalar product on qo and {tp, (p)^^ := ($ii, <&w)q is a scalar product on Qq. 
The completions which we denote by q and £2 are subspaces of L^(R^) and S) 
respectively. The quadratic form {tjj, Wzip) naturally extend to £2 and yields the 
self-adjoint Hamiltonian Bz- 

For completenes we note that for Z < Zc Equation (l5lT) shows that the norm || • || q 
is equivalent to the "Sobolev" type norm ||m||vKo ■= n(V" + f7i+ + l)|(0m^„, -u)!^. 

To prove the claimed unboundedness we pick a family of trial sequences a de- 
pending on an integer N € N - for readability we refrain from indicating this 
explicitly - given by 



(33) 



if n > TV. 



We compute the expectation of the two summands v"^ , a £ {0, 1}, of the potential 
energy and obtain 



AT /N-k+1 



(34) 



N 



-Yik + a + l) J r , +0(7VO) 



= (! + f log(A^) + 0(7V") 

for large N . Thus 

(35) £o[a] = (a, ta) ~ Z{a, ^^-±^a) = 2(1 - Z/Z^) log(iV) + 0(7V"), 

i.e., the form is unbounded from below for Z > Zc- □ 
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Appendix A. Useful Facts on the Weyl Operator with Homogeneous 

Magnetic Field 

For the convenience of the reader and for fixing the notation we collect in this 
appendix some facts related to the Weyl operator Wq with homogeneous magnetic 
field. 

We write (x) for the n-th generalized Laguerre polynomial with parameter a 
(Hochstrasser [121 Formula 22.2.12]). For m £ Z and n E Nq this allows to define 
the functions 



(36) (f>Tn,n{z) 



e 2 



k^^z^"^\Ll^\(zz) ifm>0 



7r(n+ |m|)! [e-^^^zl^lLk"' (zz) if rn < 0. 
In polar coordinates z = rexjp{iip) these functions are written as 



(37) ^^Ar,^) = ^^(;7:p|^^-^V'"'4"'(^^)e'- 

where - in abuse of notation - we use the same notation despite the change of 
coordinates. Note that these functions form an orthonormal basis of L^(M^) which 
follows from the fact the (27r)~^/^ exp(im(/?) are an orthonormal basis of i^(0, 27r) 
and for every fixed to G No the generalized Laguerre polynomials L™, n G No, 
under suitable renormalization, are an orthonormal basis of L^((0, cxd), r™e^''dr) 
(Hewitt [n]). 

Using the recursion relations [TH 22.7.29-32] of the generalized Laguerre poly- 
nomials i™ and L™'(x) = —L^^l which is immediate from the definition, we have 



(38) d*(j),n,n = 2isgn(m)^n + TO+ + l(l),n+l,n+e{-rn) 

(39) d(j)jn,n = -2isgn(TO - l)y^n + m+(l)„,_i^n-e(-(m-l)) 

(40) dd*4)„i,n = 4(n -I- TO+ -|- l)(f>m,n 

(41) d*M*rVm,« = isgn(TO)0„+i,„+e(_,„) 

for n e N, TO G Z, where - as usual - to+ max{0, to}, 

Jl ifx>0 jl ifx>0 

'(")^=\o if.<o""^'^"(")^=\-i if.<0 

(see also [20l Section 7.1.3]). Note that this solution is related to the non-relativistic 
Schrodinger equation with homogeneous magnetic field in two dimensions (Fock [7]). 
The angular momentum operator L is given as 

(42) L := xip2 - X2P1 = zd - zB = jd^. 

Writing the 0m, n in spherical coordinates easily shows that they are eigenfunctions 
of L with eigenvalue to, i.e., 

(43) i0m,n = 'rn(j)m,n- 

Since the (j)m.n form an orthonormal basis, the eigenvalue equation (j45p implies 
that \d*\ is invertible and is an isometric operator. The spinors 

for n G No and m G Z form an orthonormal basis of as shown in the proof of 
Lemma [TJ 

Using dSHl) through (gl]) we find 



(45) WoV'm,n = 2^71 + TO+ + l-0™,„, n G N, TO G Z, 



10 



T. MAIER AND H. SIEDENTOP 



i.e, for fixed m G Z, the spinor ?/'m,n is the ?i-th eigenvector of Wq on the positive 
spectral subspace F). 

Tlie total angular momentum operator J on L^(R^) is given as 

(46) ' = ^+1(1 -1 
The Formulae (gl]), (gS]), and (gH) imply 

(47) JV'm.n = (m + l/2)V'm,n- 

In fact, 

(48) [Wq, J] = and [wq, L] = 0. 

Eventually note, that $*J$ = L + 1/2 where $ is the unitary map defined 
Lemma [TJ This is the reason why it is equally natural to label the basis by the 
orbital angular momentum quantum number m as to label it by the total angular 
momentum quantum. We choose m since the formulae are easier to handle. 

Appendix B. Some Useful Facts on Gamma and Related Functions 
The Gamma function 

T(z) / z*e~* — 
^ ' Jo t 

is obviously positive on the positive half axis where it is also analytic and 
log-convex (see, e.g., Rudin [TBI Theorem 8.18]). 

A useful combination of Gamma functions is the Pochhammer symbol 

<«) «• - ^ 

which is a meromorphic function in both variables z and a. 

Lemma 4 (Gautschi ^10, Formulae 6 and 7]). For x E M-(_ and < s < 1 

(50) (. + ir.<I<f±il<..-, 

Note that Gautschi claims the inequalities for a; G N only. However, his proof is 
valid also for x G R4.. 

Furthermore, we note the reflection formula [H Formula 6.1.17] which states that 
for < 5Rz < 1 

(51) r(z)r(i-z) = -^. 

sin(7rzj 
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